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Abstract —In this paper, we propose a feedback reduc¬ 
tion scheme for full-duplex relay-aided multiuser networks. 
The proposed scheme permits the base station (BS) to 
obtain channel state information (CSI) from a subset 
of strong users under substantially reduced feedback 
overhead. More specifically, we cast the problem of user 
identification and CSI estimation as a block sparse signal 
recovery problem in compressive sensing (CS). Using 
existing CS block recovery algorithms, we first obtain 
the identity of the strong users and then estimate their 
CSI using the best linear unbiased estimator (BLUE). 
To minimize the effect of noise on the estimated CSI, 
we introduce a back-off strategy that optimally backs- 
off on the noisy estimated CSI and derive the error 
covariance matrix of the post-detection noise. In addition 
to this, we provide exact closed form expressions for 
the average maximum equivalent SNR at the destination 
user. Numerical results show that the proposed algorithm 
drastically reduces the feedback air-time and achieves a 
rate close to that obtained by scheduling schemes that 
require dedicated error-free feedback from all the network 
users. 

Index Terms —Feedback, Decode-and-Forward, Full- 
Duplex Relaying, Compressive Sensing. 

I. Introduction 

In a relay aided networks, a relay station (usually low 
powered) is deployed to assist communication between 
the base station (BS) and a destination user. Relaying 
techniques can be classified, based on their forwarding 
strategy and required processing at the relay terminal, as 
decode and forward (DF) or amplify and forward (AF) 
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|[5l . Relay transmission for both AF and DF relaying 
can be performed in half-duplex (HD) (e.g. 0) ox full- 
duplex (ED) (e.g. Q ) mode. In HD mode, the BS and 
the relay transmit on orthogonal channels, whereas in 
FD mode, the BS and relay share a common channel and 
the relay transmits and receives simultaneously over the 
same channel. Thus, half-duplex relay schemes induce 
a throughput loss due to the pre-log factor 1/2 lH. Full 
duplex relay schemes are shown to improve the network 
throughout by eliminating this pre-log factor Q. 

In a network with large number of users, the BS can 
exploit the multiuser diversity to maximize the network 
capacity. This is achieved by requesting all users to 
feed back their equivalent SNR (BS-relay-user) to the 
base station (BS). Clearly, this is expensive in terms 
of spectrum utilization and results in a large feedback 
overhead, especially when the number of users is high 
However, channel state information (CSI) mis¬ 
match that could arise from both channel estimation 
errors and quantized feedback results in a loss in the 
achievable rate. Therefore, for practical scenarios, feed¬ 
back strategies should be able to reduce the overhead 
while satisfying certain performance guarantees. In fact, 
to reduce the feedback, it is possible to request only a 
few strong users to encode their SNR on non-orthogonal 
codewords. This creates a sparse user regime, where 
sparse approximation algorithms (see e.g. ifTTl ) are then 
applied to recover the identity and SNR of the strong 
users. The work in ifT^ uses Compressed Sensing (CS) 
to reduce the feedback load of single-input-single-output 
(SISO) networks, while other work e.g. ifTTl - lfT^ apply 
CS to reduce the feedback overhead for MIMO networks. 
In this paper, we apply CS to reduce the feedback 
overhead of a relay-aided multiuser network. We exploit 
the feedback nature of FD relays to pose the feedback 
problem as a block CS recovery problem. Block CS 
recovery leverages the a priori information of the signal 
block size to better differentiate true signal information 
from recovery artifacts. This leads to a more robust 
recovery ll22l . To the best of our knowledge, this is the 
first paper that proposes a CS-based feedback algorithm 
for multi-user FD relay-aided networks. We focus on 
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networks that are comprised of a BS serving users in 
their respective coverage areas. It is assumed that there 
is no direct link between the BS and the users, and hence, 
communication must take place via a relay terminal. 
Examples of such a topology include cell-edge users, 
users that are shadowed from the BS, e.g. users in a 
shopping complex or an air-port, or users covered hy 
a pico or a femto cell as the case in heterogeneous 
networks ll20l 11211 . To exploit the multi-user diversity, the 
BS requires feedback of CSI from all the users. Based on 
the feedback the BS receives (via the relay), it schedules 
the user with the best BS-relay-user channel condition. 
This feedback requirement poses two immediate chal¬ 
lenges. Firstly, CSI feedback (users-relay-BS) generates 
a great deal of feedback overhead (air-time) that could 
result in significant performance hits if the number of 
users becomes too high. Secondly, the fed back channel 
information is usually corrupted by additive noise. In this 
paper we addresse the above challenges and provide the 
following contributions 

• We propose a CS-based full-duplex feedback algo¬ 
rithm that reduces the feedback overhead of multi¬ 
user relay aided networks. To the best of our 
knowledge, this is the first paper that introduces 
full-duplex feedback and accounts for the relay self¬ 
interference on the feedback link. 

• We account for the feedback noise and refine fhe 
esfimafed equivalenf SNR, and optimally backing 
oj^ on fhe noisy equivalenf SNR esfimafes. 

• We derive fhe posf-CS defecfion noise variance. 

• We derive fhe average equivalenf SNR of fhe besf 
user when fhe relay employs HD and FD feedback. 

The remainder of fhe paper is organized as follows. In 
Section m we infroduce fhe sysfem model. We describe 
fhe proposed feedback algorifhm in Secfion |III] and 
evaluafe ifs performance in Secfion JV] In Secfion |Vl 
we presenf some numerical resulfs prior fo concluding 
our work in Secfion |Vll 

Notations: Throughout this paper, we use the fol¬ 
lowing notations : Fx {■), {.), px{ ) and E(X) 

stand for the cumulative density function (CDF), the 
inverse CDF, the probability density function (PDF), and 
the expectation of the random variable X respectively. 
We denote by P(.A), the probability of the event A. 
Matrices are denoted by bold capital letters, rows and 
columns of the matrices are referred with lower case bold 
letters. We use the subscript c and r to refer to a matrix 
column and row respectively, for example ac,j and j 

'The term back-off denotes an SNR reduction, applied by the BS 
to the estimated SNRs, to reduce the likelihood of an SNR over¬ 
estimation error. 



are respectively the ith column and row of the matrix A. 
We denote by Aj, the submatrix of A whose columns 
are in X. tr (.) and [.]* respectively denote the trace and 
the transpose operators on matrices. To refer to a matrix 
A as a positive semi-definite matrix, we use the notation 
A ^ 0, consequently A ^ B if and only if A — B ^ 0. 
We also denote by T (.), 7 (.,.) and T (.,.), the Gamma 
function, the lower incomplete Gamma function and the 
upper incomplete Gamma function respectively. 

II. System Model 

In this section, we introduce the downlink and feed¬ 
back (uplink) transmission models for the relay-aided 
multi-user network. 


A. Downlink Transmission Model 


We consider a network with N cell-edge users served 
by one BS via a single FD relay (no direct path exits 
between the BS and destination users). As shown in 
Fig. m both the BS and the users are equipped with 
a single antenna, whereas the relay is equipped with 
one transmission antenna and one receive antenna. The 
two relay antennas operate simultaneously over the same 
frequency, i.e. reception and transmission occur at the 
same time over the same frequency. We assume that a 
direct link exits between the BS and the relay, hence, 
the channel from the BS to the relay / is modeled by 
a Rician distribution with fading parameter K = 
where 2 c 7 ^ is the average power in the Non-Fine Of 
Sight (NFOS) multipath components and 6 ^ is the power 
in the FOS component HJ. For simplicity, the BS-relay 
channel can be approximated to follow a Nakagami 
fading channel parameterized hy P = iP' + and 
fading parameter d = Therefore, 


vr 


T = 


7 {d, f) 

T{d) 


( 1 ) 


where 6 = ^. The channels from the relay to the N 
users, {gn)n=i 2 n assumed to be independent and 








3 


Feedback period {LTms) 


Transmission period 


E 


Coherence Time (Tc) 


Figure 2. Structure of a frame of duration Tc composed of feedback 
period and transmission period. 


identically distributed (i.i.d.) zero-mean complex Gaus¬ 
sian random variables with variance < 7 ^. To maximize 
the downlink link rate, the destination user is selected 
based on the following rule 

n* = argmax 7 „, ( 2 ) 

n 

P \q P 

where 7 ^ = is the instantaneous receive SNR at 

the reth user which is an exponential random variable 

PrCr^ 

with mean 7 = and Nq is the additive noise 

variance. The relay may apply cancellation techniques 
||2^ . ||25]| to (partly) eliminate this self-interference. Usu¬ 
ally this involves estimation of the loop channel h 
and then subtraction of the interfering signal. However, 
in practice, complete interference cancellation is not 
achieved due to non-ideal channel estimation and signal 
processing. In this paper, we model the effective gain of 
the loop channel due to the residual interference after 
cancellation as a multiplicative factor n that remains 
constant during the transmission intervaloLet (f) and 
Sn* {t) represent the original and the decoded symbol by 
the relay respectively, then after selecting the destination 
user n*, the BS broadcasts the symbol Sn*(f) which is 
received by the relay as 

yr{t) = ^/Ps/Sn-it) -h py/^Sn-{t “ 1) + Zr{t), (3) 

where Ps and Pr represent the BS and the relay trans¬ 
mission powers respectively, and the term Zr represent 
the zero mean variance Nq additive Gaussian noise at 
the relay. During the second hop, the relay will decode 
and forwar41 the symbol, and the n*th user receives 

2/n* (f) — y/'frfl'n* 'Sn* (f) T (4) 

where Zn*{t) is the zero mean variance Nq additive 
Gaussian noise at the n*th user. The transmitted symbols 
are assumed to have the same constant normalized 
transmitted power over time, i.e. E[|s(f)p] = 1. 


^This assumption can be justified by the fact that relay’s antennas 
are fixed during transmission. 

^Throughout the paper, we only consider decode and forward as 
the downlink protocol. 


B. Feedback Model 


We assume a reciprocal TDD setup so that the BS- 
relay and the relay-user channels on the uplink are the 
same as those of the downlink. As shown in Figure 2, 
the feedback channel is assumed to be slotted, shared, 
and opportunistically accessed by the users. There are 
L feedback mini-slots and each feedback transmission 
is received over a duration of one feedback mini-slot 
time Tms (see Figure 2). All channel coefficients, i.e. 
/, h and {gn)n=i 2 N’ assumed to be constant 
during the feedback and data transmission periods, (i.e. 
all channels are fixed during one coherence interval of 
duration Tc). Each user is assigned an ID which is a 
unique signature sequence vector of dimension M N 
for use on the slotted feedback channel. The choice of 
M is discussed in Section rTlT-CI The signature sequence 
vectors, are drawn from the columns of a real Gaussian 
matrix A S (with zero mean and variance 

i.i.d. entries), and are deterministically assigned by the 
BS to the users. 

Prior to feedback, the relay broadcasts a pilot from 
which the BS and the users estimate their channels with 
the relay. Only users having an instantaneous receive 
SNR (Please refer to eq. ©) higher than a predetermined 
threshold ^th encode (or multiply) their SNRs (relay- 
user SNR) with their signature sequence vectors and 
transmit the combination to the BS (via the relay) after 
applying proper uplink channel compensation, i.e. the 
nth user scales its transmission by — _ The remaining 
users remain silent or effectively feed back a “0”. Let 
X = [xi X 2 ... Tat]* denotes the fed back sparse 
vector, where is the fed back value of the nth user 
and it is given by 


{in ,ln> 1th 

y 0 , otherwise. 


(5) 


Also, let &r,m — 0 ,m ,2 S M , be 

the mth row of A. Then at the mth feedback mini-slot, 
the nth user scales its feedback information by am,n and 
feeds back the combination to the relay which simply 
forwards what it receives to the BS. The relay can feed 
back the received measurements using FD or HD modes. 
In the HD case, the relay receives the following at the 
mth feedback mini-slot 


y 


(m) 


^x -h 7 


(m) 


( 6 ) 


In the case of FD feedback, the BS will have to account 
for the relay self-interference which does not exist in the 
case of HD feedback. 


( 7 ) 

+ = 2,3, ...,ih(8) 
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where is the received signal at the relay at the 
mth feedhack mini-slot, oT/r™ is the self-interference 
residual at time m, and Zr is the relay noise component 
at the mth mini-slot. For the sake of simplicity, we as¬ 
sume that p is small enough such that \p\^ k, > J, 
where J is a constant, so that all factors multiplied hy 
\p\^ will he ignored for k > J. With this in mind, we 
can transform the recursive expression in Q and (7) to 
the following form 

min( J—l,m—1) 

^ p^a.r,m-k^ 

(9) 

min ([ .J ] — 1 ,m—1) 

+ ^ 

k=0 


The term ^ in the second term of (|9ll is due to the 
assumption of the convergence of the noise variance 
to zero for values of k > i.e. « 0 . 

After receiving the mth feedhack measurement, the relay 
forwards that signal to the BS. After equalizing hy the 
BS channel, the received mth measurement at the BS 
becomes 


min( J—l,m—1) 

^ p^a.r,m-k^ 
k=0 

min(|'.j] —l,m—l) 

+ J]; + 

k=0 


( 10 ) 


w 


(m) 


f ’ 


where is the noise term at the BS at the mth 

reception instant. 


III. Opportunistic User Selection and 
Feedback SNR Estimation 

In this section, we present the proposed CS-based 
feedback algorithm and outline the steps required to 
obtain the identity and the receive SNR of the strong 
users at the BS. We first demonstrate our approach in 
the simpler case of FID relaying and then extend it to 
the case of FD relaying. 


A. Feedback Threshold 

To apply CS theory, only a few strong users are 
allowed to feedback. This creates a sparse user regime. 
This is done by using a feedback threshold 0. The 
feedback threshold 7 th is optimized to meet a target 
scheduling outage probability Vo- The outage probability 
is defined as the probability that all users fail to report an 

"'Note that the threshold is the same for both relaying techniques 
HD and FD. 


SNR above 7 th. The scheduling outage probability can 
be calculated as 

Vo = P(7n < 7 th, for all n = 1,..., N) 

= [i^7(7th)]^ 

where S is the number of users that report an SNR above 
the threshold and S' = E [S'] is its mean value , Fly ( 7 th) = 
1 — exp(^l^i) is the SNR CDF at the users side. From 
CD, the feedback threshold can be calculated as 

^,^ = F-\VI/^), (12) 

where Vo is assumed to be very small. 


B. User ID Estimation 


For the BS to make a user selection decision, it has 
to estimate the fed back vector x. In what follows, we 
show that we can reliably estimate x and hence identify 
a set of strong users. More precisely, we show that in 
both feedback scenarios (HD and FD) recovering the fed 
back vector simplifies fo solving an underfermined linear 
sysfem. 

1) HD feedback: In fhe case of HD feedback, fhe 
feedback is communicated over orthogonal channels 
fhrough fhe fwo-hop nefwork via fhe relay. In ofher 
words, fhe relay listens fo fhe users’ feedback in one 
time slof and forwards what it gets in an independent 
time slot. Thus from ®, we have the following linear 
system : 

ys = Ax -F z, (13) 


where ys 
z = 




( 2 ) 

Vs 


yi^^ 


-F 


^ + 4^) 


^„(M) (M) 

j: \ 


and 

1 1 


_ f f 

Obviously, the noise vector is uncorrelated with 
covariance matrix Sg = E [zz*]. To be able to 


determine we need to evaluate E 




. We have H 






1 

id-l)9' 


(14) 


Thus, Sg = (^Nq + Im- With this in mind, we 

can recover the support of x using regular compressive 
sensing recovery algorithms such as the least absolute 
shrinkage and selection operator (LASSO) outlined in 

CD. 


^Note that for K > 0, d = 


1-F 


2if+l ’ 


SO that d > 1. 


2if+l 
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2) FD feedback: In this case, we make use of the 
structure in the received signal in (fTOl) and then ap¬ 
plies the theory of Compressive Sensing to recover x. 
Interestingly, the problem can he formulated as a block 
sparse recovery problem where the unknown vector is 
a block sparse vector with block size J. We denote by 
X(n) = [xn,pXn,-;P'^~^Xn\ , u = 1 , 2 , and by 


- 

^c,n — 




0,0, ...,0,ac,n (1 - M - j) 


Then, based on (fTOl) . we have the following linear 
system 


where B = 

( 1 )’ ( 2 )’ 

n = 1,2, ...,A^ and z = [zi Z 2 

vector at the BS, where 


ys = Bv + z, (15) 

^ [B( 1 ),B( 2 ),...,B(jv)], V = 

and B(„) = aja^-^], 


zm\ the noise 

m = 1 


Zm. — 




(16) 

It is not difficult to see that the noise at the BS is 
correlated with correlation matrix = E [zz*]. In the 
case where J = 3, we have H 


S. = 


oti 03 0 

as 02 03 

0 03 ■ ■. 
0 . 


■ • ^3 

O3 O2 


(17) 


where oi 


- ^o(l + p^), 012 

No (1 + p"^ + and 03 = pNo. 


As shown in (fTSl) . the unknown vector v is block 
sparse with block size J and sparsity S. Therefore, 
it is possible to apply the results in ll 22 l to reliably 
recover v. As shown in ll22ll . for reliable CS recovery, 
the number of measurements should scales as M = 
C (JS" + S'log JA^/S"). Basically, the BS performs CS 
recovery (e.g. the CoSaMP algorithm in |[22l ') on the 
resultant measurements vector and estimates the ID of 
the strongest user. If at least one user is detected, the BS 
performs a subsequent SNR estimation and refinements 
to be discussed below, otherwise an outage is declared. 
In the following subsection, we summarize an important 


®Our analysis is applicable for any value of J, but for presentation 
purposes, we restrict the analysis to J = 3. 


result from CS theory for recovering the sparsity pattern 
of a block-sparse vector in a noisy setting. 


C. Block-Sparse Signal Recovery 

The theory of compressed sensing permits efficient 
acquisition and reconstruction of a sparse signal (through 
multiplication by an appropriate random matrix) from 
only few measurements. For a sparse vector with size 
N and sparsity S, it has been shown that robust signal 
recovery is possible from M = 0{S\ogN/S) fV2\ . 
One would expect that additional structure imposed on 
the unknown sparse signal would substantially reduce 
M without sacrificing fhe recovery performance. One 
imposed sfrucfure could be fhe block-sparsity structure, 
where the locations of the significant entries in the 
sparse signal cluster in blocks. In the literature, the 
concept of recovering a block-sparse vector with reduced 
number of CS measurements has been studied in |[T3l 
and m. However robustness guarantees are restricted 
to either conventional CS sparse signals or recovery with 
noiseless measurements and do not have exact bounds for 
the required number of CS measurements. Recently in 
Baraniuk et al. showed that robustness guarantees 
can be achieved with M = C{JS Slog JN/S) CS 
measurements, where C is a positive constant and J 
is the cluster or the block size, which is a substantial 
improvement over M = C {JS log N/S) that would 
be required by conventional CS recovery algorithms. 
For more analytical details about block-spare recovery, 
readers are referred to ini-iEa and references therein. 


D. User SNR Estimation 

1) HD feedback: After CS recovery, the BS obtains 
information on the location of the non-zeros in x or 
equivalently the support of x denoted by S, where 
|5| = S is the cardinality of S. Therefore, the linear 
system in ([T3] ) can be rewritten as 


Ys = A 5 X 5 -f z 


(18) 


We can now apply least squares (LS) to estimate the 
entries of X 5 as follows 


x<s — (A5A5) A^Ys 
= X5 -f ols , 


(19) 


where 0^5 is the estimation error vector. Upon condi¬ 
tioning on A 5 , vector e^s is Gaussian since it results 
from a linear operation on Gaussian random variables. 
The covariance of 6^5 is given by (A^A^) ^ . As M 
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grow to infinity while S is fixed, (A^A^) ^ becomes 
close to its mean ll2^ in the sense that : 

The convergence holds in the almost sure sense and 
shows that the variance of the error is almost the same for 
all vector entries. In light of this result, we will consider 
for simplicity that e is is Gaussian vector whose entries 
are indendenent Gaussian variables with zero-mean and 
variance al^g = m^Nq. 

2 ) FD feedback: Similarly, in the case of FD feed¬ 
back, the BS recovers the support of v denoted by X, 
where |X| = JS. Therefore, the linear system in (fTSl) 
can be rewritten as 


ys = BiVi-hz, (21) 


Since the covariance matrix Sg is correlated, one can 
apply the Best Linear Unbiased Estimator (BLUE) 0 to 
estimate the entries of vj as follows ll2^ 

vx=(B^5]riBx)’'B^SrV. ^^2) 

= VI -I- eBLUE, 


where gblue is the BLUE estimation error. Similarly 
to the HD feedback scenario, upon conditioning on 
Bi, vector gblue is Gaussian with zero-mean and 
covariance (BjS- ^Bi) . When M tends to infinity 
with S and J fixed, we can prove that: 

Lemma 1 : Assume that M —)• -|-oo with S and J fixed. 
Assume that limsup^y^ 11^2^11 < +oo- Then, 

(B^SriBi) 0. (23) 

Proof: See Appendix A for proof. ■ 

Again, we note that the above convergence suggests con¬ 
sidering the noise at the output of the BLUE estimator as 
a Gaussian random vector whose entries are independent 
with zero-mean and variance ct'^blue given by: 


^e,BLUE — 


M 


tr (sri) 


( 24 ) 


E. SNR Back-Off 

As stated in ([T^ and (l 22 l) . the LS and BLUE esti¬ 
mators are noisy which means that the estimated SNR 
can be higher or lower than the actual one. This is 
problematic, since an estimated SNR higher than the 
actual one results in a transmission rate higher than the 
maximum rate the end user can support. To deal with 

’when the noise is white as in eq. idl, LS and BLUE are 
equivalent. 


this, we back-off on the estimated noisy SNRs. Erom 
(fT^ and (l22l) . each entry of and vj respectively can 
be represented in a scalar equation as 7 = 7 -I- e, where 
7 and 7 stands for the actual and the estimated SNRs 
respectively, and e is the Gaussian error. Our back-off 
strategy simply subtracts an amount A from 7. Hence, 
7 becomes 

7 = 7 -p e — A ( 25 ) 


To characterize the outage when the estimated SNR is 
higher than the actual one, we define ry as fhe back-off 
efficiency which is simply given by 


ry =P (7 < 7) 
=P (e < A) 



( 26 ) 


where cje is either the LS or the BLUE error variance. 
Obviously, the performance of the proposed algorithm 
will depend on A. A low value of A may result in an 
outage while a high value will result in a low rate. The 
optimal value will be derived in the following section. 


IV. Performance Analysis 

The performance of the proposed feedback algorithm 
will be evaluated based on three performance criteria: 
i) the feedback load, ii) the achievable rate, and iii) 
the achievable throughput. Obviously, each of these per¬ 
formance criteria will depend on the feedback relaying 
mode, i.e. HD or ED feedback. HD feedback results 
in lower noise variance at the BS at the expense of a 
larger feedback overhead, while ED feedback results in 
a higher noise variance at the BS (due to the relay self¬ 
interference) but results in a lower feedback overhead. 


A. Feedback Load 

The feedback load L is defined as the total number 
of feedback mini-slots required for the BS to make a 
user scheduling decision. In other words, L is the total 
number of measurements required to have robust CS 
recovery at the BS and it is given by 

Lhd = 2 CS\ogN/S 
in the case of HD feedbaclH and 

Lfd = C{JS + S\ogJN/S) , 

*In HD feedback (J = 1), the relay receives and forwards the 
users’ feedback information using two orthogonal channels. There¬ 
fore, the BS requires twice the number of mini-slots for feedback 
reception when compared to ED feedback. 
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in the case FD feedback |[22l . The term S represents the 
average number of users that feedback and it is given by 


N 


n=l 


N 


\N—n 


= N i-n 


l/N 


7^ = E 


where 


log (1 + 7 eg - A){1 - Vo) [1 - Q 


leq — 


mm 


Ps\f\^ 


) Tn* 


exp ( —) log (1 + 7 eg - A*) 


+ Q 




= 1 


^7e, ( 3 ;) =1 - 


= 1 - 


7 (d, t) 


1 - 


r(d) 


1 - 1 - e” 


N' 


where A = + Aq. Therefore, (x) is given by 

.5^7e, (X) 


S = Y,n(^' ]{l-F^{jth))^F,{^thy 

\ n / (27) 


Pie, (X) 


dx 

N _£ / / , AxA 


B. Achievable Rate 

The achievable rate for the HD and FD cases can be 
expressed as 


+ 


1 - 1 - e' 


Gi{x) 

N' 


r(d) \e 


d-l-^ 


— \ X e 0 


G2{x) 


(33) 


Hence, 


(28) 


(29) 


7eg = / XPey^^ (x) dx 

Jo 

roo 


e - -e,LS in the 
e BLUE Otherwise. Using 
can be upper bounded as □ 


PrP^ + No' 

where n* is defined in eq. Q, and 
case of HD feedback, 
the Jensen’s inequality, 

7^ < log (1 + nieq] - A) (1 - Vo) (^1 - Q . 

(30) 

It is clear from (l30l) that the achievable rate depends 
on the back-off value A. To obtain the optimal value of 
A, we maximize the achievable rate with respect to A. 
To do that, we differentiate the upper bound on TZ with 
respect to A and equate the result to zero. We denote by 
A*, the optimal back-off, then A* satisfies 

'l +%g- A* ^^_ ( 

A 


= / xGi [x) dx+ xG2 {x) dx 


(34) 


Pi 


P2 


To be able fo derive fhe infegrals in (1341) . we use fhe 
fact that 

and fhaf {u,ax) dx = 

Jll+pyJ 2 Fi(l,p + v-p + l-^^ EH, where 
2 U 1 (.,.;.;.) is fhe Gaussian hypergeometric function. 
Therefore, the first integral in (l34l ) can be evaluated as 

N\^ (N-l\ y (\ n + r 

2^d{d +1)9<^ y n y 7 


X 2 Ui ( 1, d -I- 2; 3; 


jn+l) /7 
X/9 + {n + 1) /7 


(31) 

where, 7 eq = E [7eq]- Then, the value of A* that satisfies 
(|3T]) is used in (I2H . 

If now remains fo derive fhe average end-fo-end SNR 
af fhe sfrongest user, i.e. E[ 7 eq]. To achieve fhis, we 
firsf need fo derive fhe PDF of ■jeq, Pie, (^)- Since 
7 eq is simply fhe minimum of fwo independenf random 
variables, ifs CDF can be derived as follows 


(35) 

Similarly, fhe second integral in (l34l ) can be calculated 
as 

-d-l 




fN' 


n 


n=l 


(- 1 ) 


n+l / ^ I ^ 
0 7 


Finally, 


leq = Ft + P2 


(36) 

(37) 


(32) 


®Note that the inequality comes also from the fact that we didn’t 
include the CS recovery probability in the rate expression. This 
probability is not available in closed forms for block-CS recovery. 


C. Achievable Throughput 

In the previous section, it was assumed that the 
amount of air-time reserved to feedback is negligible 
as compared to the transmission time. This assumption 
does not give much insight for practical scenarios. In 
this paper, we define fhe achievable fhroughpuf as fhe 
number of fransmiffed bifs per unif fime (bps/Hz). The 
nefwork fhroughpuf can be explicifly given by 

(2c “ LTfns) 

~Tc (38) 


T = n- 


= 'R{1-Lt) 
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Figure 3. Average Receive SNR versus the number of users (N), 
0-2=0 dB, cTg = 5 dB, No = -15 dB 


Figure 4. Average Rate versus the number of users (N), = 20 

dB, cr2 = 0 dB, cTg = 5 dB, and Vo = lO"^. 


where Tc is the channel coherence time, T^s is the time 
needed to transmit one feedhack mini-slot and r = 
is the normalized mini-slot time (see Figure 2). 

V. Numerical Results 

For all simulations, we set (7 = 2 and we use the 
CoSaMP algorithm proposed in mB for CS recovery. 
Unless otherwise specified, we set: Pg = Pr = I, 
62 = 20 dB, 0-2 = 0 dB, 0-2 = 5 dB, Nq = -15 
dB, Pq = 0.01, and J = 3. Throughout this section, 
we compare the proposed scheme (scheme. 1) with the 
following schemes: 

• Noiseless Dedicated Feedhack (scheme. 2): 
All users feedhack in independent noise-free, 
interference-free feedhack channels (via the relay). 

• Random User Selection (scheme. 3): For each trans¬ 
mission period, a user is randomly selected and 
allocated the downlink channel resource. 

• Proposed without hack-off (scheme. 4): The same 
as scheme. 1), hut without any hack-off (SNR 
refinement) strategy. 

In Figure [3l we compare the theoretical result obtained 
in eq. (ITTI) with the numerical (simulated) average SNR. 
For all values of the average LOS component the the¬ 
oretical and the numerical results exactly match, which 
validates the formula in eq. (iTTl) . It is also clear from 
Figure [3] that increasing b'^ increases the average receive 
SNR as well. However, we notice that the average receive 
SNR converges for b'^ > 10 dB, this can he explained 
hy the fact that the receive SNR is constrained hy the 
SNR of the second hop. In Figure IH it is clear in one 
hand that scheme. 2 offers the best performance in terms 

’°CoSaMP recovery algorithm that we use to recover v is available 
at http://dsp.rice.edu For HD feedback, the group size becomes 1. 



Figure 5. Average Feedback Load versus the number of users (N) 



Figure 6. Average Throughput versus the number of users (N), r = 
1/600 


of average rate, which is expected since perfect channel 
state information (CSI) is available at the BS. In the 
other hand, both versions of the scheme. 1 (HD and 
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Figure 7. Average Throughput versus the number of users (N), r = 


FD feedback) offer a comparable performance with a 
gap of 0.42 bps/Hz as compared to scheme. 2, where 
the HD version is slightly better than the FD version 
as it is free from interference. Figure |4] also shows that 
the performance of the proposed scheme without back-of 
(scheme. 4) suffers a rate-hit as compared to scheme. 1 
and scheme. 2, and hence, the BS has to perform a back¬ 
off strategy in order to account for the feedback noise. 
Random user selection (scheme. 3) is shown to yield the 
worst performance since no feedback is communicated 
to the BS. 

In Figure [51 we evaluate the performance in terms of 
feedback load, where we assume that the amount of time 
needed for channel estimation is negligible compared to 
the feedback air-time. As shown in Figure [51 scheme. 1 
consumes much less feedback as compared to scheme. 2 
especially at high number of users N. This is expected, 
since scheme. 2 requires the feedback to grow linearly 
with N, whereas scheme. 1 requires the feedback to grow 
logarithmically with N. The feedback load of scheme. 1 
with HD feedback is worse than that of scheme. 1 with 
FD feedback since feedback is performed over the two 
hops. 

To have more insight, we simulate the network 
throughput for r = 1/600 in Figure 6. As shown, 
scheme. 1 offers the best throughput, and scheme. 2 
yields the worst throughput since it requires a large 
amount of feedback. The throughput of the HD version 
of scheme. 1 is shown to deteriorate with the number of 
users. The reason for this is that HD feedback requires 
more feedback time as the number of users grow when 
compared to FD feedback. In Figure [71 we decrease r 
to be r = 10“'^ (i.e. higher coherence time). We first 
notice that at low number of users (N < 250), scheme. 
2 performs better than scheme. 1 (due to the reduced 


value of r) and the HD version of scheme. 1 performs 
slightly lower than scheme. 1 (FD). At high number of 
users, scheme. 1 (FD) performs the best. Again this is 
due to its lower feedback requirements. 

VI. Conclusion 

In this paper, we proposed an opportunistic CS-based 
feedback algorithm for a relay-aided multiuser network. 
Instead of allocating a feedback channel for each user, all 
users are allocated a pool of shared channels for feedback 
transmissions. We derived the end-to-end receive SNR at 
the destination user and the post-CS detection/refinement 
error variance. Moreover, we considered both HD feed¬ 
back and FD feedback and showed that FD feedback 
offers comparable rate performance when compared to 
HD feedback. In addition to this, FD feedback results 
in lower feedback load, especially at high number of 
users, and therefore, offers better throughput than HD 
feedback. 


Appendix A 
Proof of lemma 1 

Let H = Matrix H is of finite size JS. It 

thus suffices fo study the convergence of the elements of 
H. Using results on the convergence of quadratic forms 
ll27l Lemma5,Lemma4], we can establish that the off- 
diagonal entries of H given by 

(ay*5]riai,fc,V(n,j)/(fc,0 


converge almost surely to zero, while the diagonal entries 
converge almost surely to their mean, i.e. 


(ain) 


— r! ■ 


M^+cx) 


A 0. 


with dj^n = IE 


*-c,n 


^^i,n reading as: 


In the expression of dj^n, the dominant term is 
■^tr while the second term converges almost 

surely to zero, since 


1 

M 



M 


< ;^hmsup||5]- I 


We therefore obtain: 

H„i _ Mljs 

tr (Sl- j M^-\-oo 

In Figure [8j we plot the exact and almost sure limit 
(lemma 1) of the noise variance versus M/S. From the 
figure we notice that the lower bound on is 

tight for moderate to high ratio of M/S. This concludes 
the proof of the lemma. 
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Figure 8. BLUE error variance: Exact versus theoretical limit 
(Lemma 1). 
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